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We investigate several factors controlling the physics of hybrid structures involving ferromagnetic 
/**\ , domain walls (DWs) and superconducting (S) metals. We discuss the role of non collinear mag- 

netizations in S/DW junctions in a spin g> Nambu ® Keldysh formalism. We discuss transport in 
S/DW/N and S/DW/S junctions in the presence of inelastic scattering in the domain wall. In this 
s ! , case transport properties are similar for the S/DW/S and S/DW/N junctions and are controlled by 

sequential tunneling of spatially separated Cooper pairs across the domain wall. In the absence of 
inelastic scattering we find that a Josephson current circulates only if the size of the ferromagnetic 
region is smaller than the elastic mean free path meaning that the Josephson effect associated to 
£j . crossed Andreev reflection cannot be observed under usual experimental conditions. Nevertheless 

a finite dc current can circulate across the S/DW/S junction due to crossed Andreev reflection 
associated to sequential tunneling. 
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I. INTRODUCTION 



A simple way of obtaining correlated pairs of electrons in solid state devices is to extract Cooper pairs from a BCS 
superconductor. Devices based on this principle have focussed an important interest recently. For instance entangled 
' pairs of electrons can be manipulated in double dot experiments^. Other devices involving a larger number of quantum 
dots have been proposed recently as a quantum teleportation experiment 2 . Devices involving several ferromagnetic 
electrodes connected to a superconductor have been investigated recent!}* 3 * 4 *^. Noise correlations can also provide 
useful information about quantum entanglement*. 
. Many phenomena are involved in the proximity effect at ferromagnet / superconductor (F/S) interfaces. For 
(*C) ' instance it is well established that the pair amplitude induced in a ferromagnetic metal oscillates in space. An 
CNJ , interesting consequence is the possibility of fabricating S/F/S 7r-junctions in which the Josephson relation is 
/ = I c sin (ip + TT^gaLifliiiiiSii&iiii^. I n F/S/F trilayers the superconducting transition temperature is larger in the 
antiferromagnetic alignment of the ferromagnetic electrodesi&ii because a finite exchange field is induced in the 
superconductor in the ferromagnetic alignment. On the other hand there exist "non local" superconducting correla- 
tions coupling the two ferromagnetic electrodes that favor Ap > Aaf (the zero-temperature superconducting order 
parameter is larger in the ferromagnetic alignment It is also well known that the superconducting transition 
temperature of F/S multilayers oscillates as the thickness of the ferromagnetic layers is increased22i2ii22i22i24i2£. Sev- 
""j" 1 i eral recent works have investigated new phenomena taking place in diffusive F/S heterostructures2Si2L22i22^2i2i. Other 
recent works were devoted to understand the interplay between Andreev reflection and spin polarization at a single 
S3 i F/S interfaco 32 ! 33 . 34 . 

In a recent article M. Giroud et &1. have proposed on the basis of experiments that the proximity effect at F/S 
interfaces could be strongly modified by the presence of Cooper pair-like states propagating along domain walls 
(DWs) 3 ^. These Cooper pair-like states correspond to pair states in which the spin-up and spin-down electrons 
propagate in a neighboring spin-up and spin-down magnetic domain. This proximity effect is not strictly speaking 
equivalent to the proximity effect a N/S interfaces. The reason is that the pair correlations induced in the N side of a 
N/S interface have entangled orbital and spin degrees of freedom^. By contrast for half-metal ferromagnets the wave 
function associated to the propagation of superconducting correlations along domain walls is given by the product 
state |e,a, |) <£> |e,/3, |), where a and represent two points in neighboring magnetic domains. Another difference 
between a N/S interface and a multiterminal hybrid structure is that the incoming electron and the Andreev reflected 
hole propagate in different electrodes in multiterminal structures. As a consequence the Andreev reflected hole cannot 
follow the same trajectory as the incoming electron. This has important consequences regarding disorder averaging. 

The purpose of our article is to investigate theoretically the mechanisms by which the Cooper pair-like state 
\e, a, |) \e, 0, J.) can propagate along a ferromagnetic domain wall and to investigate several new situations that may 
be the object of experiments in the future. In section ITU we discuss the perturbative transport formula of a S/DW 
junction in which the domain wall consists of many independent channels in parallel having a rotating magnetization. 
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FIG. 1: The device involving crossed Andreev reflection and elastic cotunneling with non collinear magnetizations. Electrode 
ending at site "a" is a ferromagnet with a magnetization pointing in the direction 9 a - Electrode ending at site "b" is a 
ferromagnet with a magnetization pointing in the direction 9t- 



To discuss this model we use the spin ® Nambu ® Keldysh formalism described in section^ For the sake of obtaining 
analytical results we restrict the discussion to the transport formula obtained within lowest order perturbation theory. 

If propagation in the ferromagnet is phase coherent then the pair state |e, a, f) <S> |e, /3, |) injected at one end of the 
domain wall can propagate to the other end. On the other hand if the phase coherence length 1^ is small compared 
to the size of the ferromagnetic region then inelastic scattering processes are strong and there are just a spin-up 
and a spin-down electron propagating independently in the spin-up and spin-down magnetic domains. There is no 
Josephson current but there exists crossed Andreev reflection taking place locally at each F/S interface, so that the 
conductance is larger in the presence of the domain wall. 

In section llVl we discuss the S/DW/N and S/DW/S junctions in a regime where transport properties are dominated 
by inelastic scattering in the domain wall. The domain wall is represented by two channels in parallel, with an opposite 
magnetization. This schematic model of domain wall is expected to capture the essential physics, and can be a useful 
comparison for more realistic studies involving numerical simulations that we plan to carry out in the future. We show 
that within lowest order perturbation the transport properties are governed by processes taking place locally at each 
interface once the summation over the different conduction channels has been carried out. The chemical potentials 
in the domain wall are determined by evaluating the current circulating through each interface and imposing current 
conservation. 

In section we consider the other situation where inelastic scattering within the domain wall can be neglected. 
In this situation a finite average Josephson current can circulate between the two superconductors of the S/DW/S 
junction only if the size of the ferromagnetic region is smaller than the elastic mean free path. This condition is 
not realized with usual ferromagnets and we come to the conclusion that there is no Josephson current under usual 
experimental conditions. Final remarks are given in section IVT1 

To end-up the introductory section we note that the theory of inhomogeneous ferromagnets with non colinear 
magnetizations in contact with a superconductor was already elaborated in Refs^*^ in connection with the long- 
range proximity effect associated to the triplet component of the superconducting condensate. In our article the 
emphasis is put on other aspects of this problem (the transport of spatially separated Cooper pairs). Both effects 
may play a relevant role in experiments. Finally a recent preprint^ appeared in which the conductance of a S/DW 
junction was calculated independtly from our work. 

II. PRELIMINARIES 

A. Spin ® Nambu ® Keldysh formalism 

The direction of the magnetization is rotating in a ferromagnetic domain wall. To describe superconducting cor- 
relations in the presence of non collinear magnetizations we use a spin ® Nambu ® Keldysh formalism22^2iiSi. The 
advanced Green's function is a 4 x 4 matrix: 

(|ct T (t'),c i>T (t)J) ({c u (t'),c tA (t)}) dc+^),Ci^t)\) ({c,-, T (f),c ilT (t)}) " 
(|ct T (i'),c+ (t)j) ({c,-,i(tV+(t)}) (jc+ (f),c+(t)j) ({ Cj -, T (t'),c+(t)}) 

<K T (*'),c u (4> ({cj,iV),ci,im ({cun*,i(t)}) <Kt(*'),cuW}> 

<{c+.(i'),c+(t)}) ({c,u(f),c+ r (*)}) <{c+ (f), C + }) ({c,-, T (f),c+.(t)}) 



Gt tj {t,t') = -ie{t-t') 
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The Dyson equation relates the Green's functions of the connected system to the Green's functions of the disconnected 
system. In a compact notation the Dyson equation takes the form G — g + g £§> £ (g> G, where the symbol (g> includes a 
summation over the sites of the network and a convolution over time variables. Since we consider stationary transport 
the convolution over time variables becomes a simple product after a Fourier transform is carried out. The Dyson 
equation for the Keldysh Green's function G + '~ is given by 



G + '~ = 



I + G 1 



T,®G J 



(2) 



where the self-energy £ contains all the couplings present in the tunnel Hamiltonian. The tunnel Hamiltonian 
corresponding to Fig. takes the form 



w = 



(3) 



The current through the link a - a is given by 



e 

2h 



where the matrix a z is given by 



Tr 



10 
0-100 
10 
-1 



I du, 



(4) 



(5) 



and the Nambu representation of the hopping matrix elements is given by t a ,, 

tp.b = t/3,b&z- 



B. Green's function of a ferromagnetic metal 



Now we give the expressions of the Green's functions of a ferromagnetic metal. We first suppose that the spin 
quantization axis is parallel to the direction of the magnetization. The Green's function takes the form 

g 1A (R,ui) 

9i,2{R,u) 

g^{R,uj) 





g(R,tu) = 

The four diagonal elements are given by 
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(10) 



where we have introduced a Fermi wave vector mismatch as well as a mismatch between the spin-up and spin-down 
Fermi velocities. The parameter ao is equal to the distance between neighboring sites on the cubic lattice. For 
generality we introduced a different mass for the spin-up and spin-down electrons, meaning that the spin-up density 
of states is different from the spin-down density of states. The local propagators are defined by 
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We also introduced phenomenologically in J7J - (|10fl an exponential decay of the correlations due to the presence 
of a finite coherence length in the ferromagnet. 1^ is usually smaller than the dimension of the ferromagnetic 
metal. In this case ferromagnetism can be treated semi-classically like in the theoretical description of the giant 
magnetoresistanca^i^Siii However Aharonov-Bohm oscillations in a ferromagnetic nanoring have been reported 
recentlji^. The inner diameter of the Fe-Ni nanoring in Refi^i is 420A and the outer diameter is 500A 

We will use in section lTTTl the expression of the local Green's functions of a ferromagnetic metal with the quantization 
axis not parallel to the magnetization. We suppose that the direction of the exchange field is rotated by an angle 9 
around the x axis. We do not incorporate a rotation of angle <p around the z axis since this rotation just introduces 
simple phase factors. The local Green's function of the rotated ferromagnet takes the form 



9lo 



ITTp 



l+Pcos9 

1 - P cos ( 
iP sm6 

iP sin 9 



-iP sin 9 

-iP sin 9 

1 - P cos 9 

1 + P cos ( 



(13) 



where p — (pj + pi)/2 is the average density of states at the Fermi level and P — (pj — Pi)/{p] + Pi) is the spin 
polarization at the Fermi level. 

We will also use in section IV Bl the expression of the full propagator g(R,oj) of a rotated ferromagnet. The Green's 
function takes the form 



g(R,cj) = 



9i,i 3i,3 

32,2 52,4 

53,1 <? 3: 3 

.94,2 g 4j4 



where the diagonal elements are given by 



9i,i = 


1, 

2(31,1 


+ 33,3) H 


- \ cos 9 (.91,1 


- 33,3) 


92,2 = 


\(92,2 


+ 34,4) H 


- - cos 9{g 2 ,2 


- 34,4) 


93,3 = 


2(33,3 


+ 3i,i) H 


- \ COS 0(53,3 


- 31,1) 


34,4 = 


1/ 

2(54,4 


+ 32,2) H 


- icos6>(.94, 4 


- 32,2), 



(14) 

(15) 
(16) 
(17) 
(18) 



where gi i, 32,2, 33,3 and 54,4 are given by (|7|l- i|lU[l . The extra-diagonal elements are given by 



31.3 = -33,i = 2 sin %3,3 - 3i,i) 

1 

32.4 = -34,2 = ^ Sm6, (32,2 ~ 34,4) 



(19) 
(20) 



C. 4x4 Green's functions of a superconductor 



The Green's function of a superconductor takes the form 

g A > R (R,u;) = 



g(R,iu) f(R,io) 

f {R,u>) g'(R,u) 

g{R,cu) -f(R,u) 

-f(R,u) g'(R,u) 



The matrix elements of the Green's function are given by 
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(21) 

(22) 
(23) 
(24) 



where we supposed that to < A. The coherence length is given by £(w) = hvF / V A 2 — uj 2 . 
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D. 4x4 Green's functions of a superconductor in a uniform magnetic field 



A uniform magnetic field h$ can penetrate in a superconductor if the superconductor is in a thin film geometry^ 
and the magnetic field is parallel to the direction of the superconducting film. The effect of the magnetic field is 
a Zeeman splitting of the spin-up and spin-down quasiparticle bands. Let us suppose that the quantization axis is 
parallel to the orientation of the magnetic field. The 4x4 Green's function takes the form 



A.R 



g+(R,w) f+{R,oj) 

f+(R,u) g' + (R, cj) 










g-{R,u) 
-f-(R,u) 






-f-(R^) 

g'_(R,w) 



(25) 



with g+(R, u>) = g{R,uj 
g'_(R,u;)=g'(R,u;-h s ), 



f h s ), g'+(R,u) = g'{R,uj + h s ), f+(R,u>) = f(R,co + h s ), g-{R,u) = g(R,u> - h s ) 
f-(R,u) = f(R,u-h s ). 



III. CROSSED ANDREEV REFLECTION AND ELASTIC COTUNNELING WITH NON COLLINEAR 

MAGNETIZATIONS 



A. Transport formula 



In this section we evaluate the transport formula corresponding to the device on Fig. ^ in which the magnetization 
of electrode "a" ("b") makes an angle 9 a (9b) with the z-axis. Using the formalism described in section [H] we obtain 
the current per conduction channel through electrode "a" to lowest order in t 0iQ and tb,p: 

I a ,a = I J dw8nX~pl(l - Pa)fL(") - eV a ) - n F (u> + eV a )] (26) 

+ | J duj4ir 2 tlt 2 p p a p b [1 + P a P b cos (9 a - 6 b )] {(g a ,^, a )) 

x — eV a ) — iif{ijJ + eVa) — «f(w - eV&) + rif(w + eV b )] 

+ ~ J dwAltHltlpaPb [1 - P a P b COS (9 a - Ob)] «/ a ,/j//J,a)> 

x [uf{^ — eV a ) — + eV a ) + np(uj — eV&) — np (uj + eVb)] ■ 

We have supposed that electrodes "a" and "b" are made of a large number of independent conduction channels in 
parallel so that we make an averaging over the microscopic phases in the propagators. Now if we consider that the 
same voltage is applied on both electrodes the conductance is given by local Andreev reflection and crossed Andreev 
reflection: 

G = 327r 2 ^t^(l - P 2 )f 2 oc + 32tt 2 ^-t 2 J 2 p a p b [1 - P a Pb cos (9 a - 9 b )} «/„, /3 / / j,a)). (27) 



B. Conductance associated to a domain wall 



Let us now consider the situation on Fig. [21 representing a S/DW contact between a superconductor and a mag- 
netic domain wall. We suppose that the ferromagnetic metal is made of a collection of independent channels. The 
magnetization is rotating inside the domain wall meaning that the angle 9 is a function of z: 9 = 9(z). We want to 
evaluate the difference G^ DW ^ — G^ between the conductances G^ DW ^ in the presence of the domain wall and G^ 
in the absence of the domain wall. To obtain the conductance we sum the contributions of the different channels (see 
Fig. EJ and we obtain 

G (DW) _ G (0) = X ^4-2^2 ( F fa D) (28 ) 

h a \ h J 

with 



f&,d)= i ^ [*± sin2 ( *(*) - °(*)) exp _ v (**y + (* - *£] , (29) 

ao J a-o J a Q (Ayy + (z a - z b ) 2 \ 2 J \ £ / 
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FIG. 2: The device involving a S/DW junction between a superconductor and a ferromagnetic domain wall. In the ferromagnet 
the local magnetization makes an angle 8(z) with the z-axis. The x- and z-axis are in the plane of the figure. The y axis is 
perpendicular to the plane of the figure. 




FIG. 3: Variation of F(£o,D) as a function of £o for different values of D. We used the domain wall profile given by (1: >( )l in a 
three-dimensional geometry. The distance between neighboring channels is ao = 1. We introduced a cut-off —L y /2 < Ay < L y , 
—L z /2 < z a ,zt < L z /2 in the expression of F(£o,D) (see Eq. 129H ). The solid lines correspond to L y = L z = 500 and the 
dashed line (almost superimposed with the solid line) correspond to L y = L z = 600. 



where £ = hvp/A is the BCS coherence length at zero energy and D is the width of the domain wall. L y is equal 
to the dimension of the junction in the y direction and we used the notation Ay = y a — yb- To obtain lj29|) we have 
supposed that the width of the domain wall is much larger than the Fermi wave-length so that we can average over 
the microscopic phase variables in the propagator f a ,p (see Eq. JUJ). 

Crossed Andreev reflection cannot take place between the channels separated by a distance much smaller than the 
width D of the domain wall because such channels have an almost parallel magnetization. Crossed Andreev reflection 
cannot take place either between channels separated by a distance much larger than the superconducting coherence 
length because of the exponential decay of the propagator f a ,p- As a consequence the value of G {DW) - G (0) is the 
largest if the width of the domain wall is small compared to the BCS coherence length. This is illustrated on Fig. [3] 
where we have represented the variation of the conductance as a function of £o for different values of D and for the 
domain wall profile given by 



9(z) = arctan 




C. Exchange field in the superconductor due to the proximity effect 

Now we come back to a system in which two ferromagnetic electrodes are connected to a superconductor. An 
exchange field can be generated in the superconductor because of the proximity effect. This was first observed in 
Ref^ in the case of insulating ferromagnets. An exchange field in a superconductor is a pair breaking perturbation. As 
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a consequence in the F/S/F trilayer with insulating ferromagnets the order parameter is larger in the antiferromagnetic 
alignment of the ferromagnetic electrodes^. This was well verified in experiments with insulating ferromagnets^il. 
The same effect is present with metallic ferromagnetsi£^ but in this case there exists also pair correlations induced 
in the ferromagnetic electrodesi&iS that can modify the value of the self consistent order parameter. 

We suppose that the magnetizations in electrodes "a" and "b" make an angle 9 a and 9 b and that an exchange field 
hs is induced in the superconductor. Without loss of generality we suppose that the direction of the exchange field in 
the superconductor is parallel to the quantization axis. In terms of the g + , f + , g_ and /_ introduced in section Til Dl 
the transport formula is found to be 



• J dLulirhipl [(fl + fl) (1 - P 2 a cos 2 9 a ) - 2/ + /_P 2 sin 2 9 a ] [n F (uj - eV a ) - n F (uj + eV a )] 



- I du>2i 
h 



P a P b COS 9 a COS 9 b 



+ 2 ((.9a,/35 / g, Q >>^ , a-Pf) sm ^aSin6' & | x [n F (uj - eV a ) - n F (u + eV a ) - n F (u> - eV b ) + n F (uj + eV b )} 
nXtjpaPb { [((ft f jfl a )) + (if-pfpj)] [1 - PaPb cos 9 a COS 9 b ] 
((fi.,pfl3, a )} p a p bsm9 a sm9 b j x [n F (u> - eV a ) - n F (uj + eV a ) + n F (uj - eV b ) - n F (uo + eV b )] 



- I dw2 



(31) 
(32) 

(33) 



The term (|31|l corresponds to local Andreev reflection. The term l|32|l corresponds to elastic cotunncling and the 
term (|3^|) corresponds to crossed Andreev reflection. The term pgt a )) corresponds to a process in which a 
spin-up electron travels from electrode a to electrode b and comes back to electrode a as a spin-up electron. The term 
((da p9p a)) corresponds to a process in which a spin-up electron travels from electrode a to electrode b, undergoes a 
spin precession in electrode b and comes back as a spin-down electron traveling from electrode b to electrode a. 

Replacing the propagators involved in Eqs.JjHJ, (|32() and (|33|l by their expressions given in section Hi CI leads to the 
transport formula to lowest order in hs and w. 
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+ -jp2 P a Pb cos (9 a + 9 b ) \ x [n F (uj - eV a ) - n F (oj + eV a ) - n F (uj - eV b ) + n F (oj + eV b )] 



+ 2TT 2 t 2 a t}p a p b 



\2irR, 
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a,f3 



exp 



2R 



a,/3 
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[1 - P a P b cos {9 a - 9 b )} 



(34) 
(35) 

(36) 



x [n F (uj — eV a ) — n F {uj + eV a ) + n F (u> — eV b ) — n F (uj + eV b )] 



where k F Rq — ir/2 is the ultraviolet cut-off used to define the local propagator involved in local Andreev reflection. 
We see that the crossed Andreev reflection term given by (|35() is not identical to the elastic cotunncling term given 
by H3tj|) . This shows that the symmetry between elastic cotunneling and crossed Andreev reflection is broken by the 
exchange field in the superconductor. This can be illustrated by considering that electrode "b" is a normal metal: 
P b = 0. The crossed conductance at zero voltage is finite if the exchange field hs in the superconductor takes a finite 
value: 



81 

(*a,b = q^- = 47T t a tpp a pb 



«0 



2nR a ,j3 



exp 



2Rg.J3 



hs 
A 



(37) 



By comparison we have G a , b = if hs = because of a cancellation between the crossed Andreev reflection and 
elastic cotunneling conductances. We thus see that a crossed Andreev reflection experiment with a ferromagnetic and 
a normal metal electrode can give information about the existence of an induced exchange field in the superconductor. 
We see also from Eqs. (|31fl . (|32() and (|33|l that there is no precession of the electron spin around the direction of the 
exchange field in the superconductor. The absence of spin precession in the superconducting case can be contrasted 
with the metallic case (see Appendix lA")l . 
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FIG. 4: The device considered in section llVl In (a) there is no domain wall in the junction. In (b) a domain wall is pinned in 
the junction. 
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FIG. 5: The device considered in section ft V Ul in which two single channel electrodes representing two magnetic domains are 
inserted in between a superconductor and a normal metal. The two electrodes ending at sites a and a 1 and sites b and b' are 
ferromagnetic. In sections II V PI and W\ we suppose that the electrode containing the sites a' and is superconducting. 



IV. SEQUENTIAL TUNNELING OF COOPER PAIRS THROUGH A MAGNETIC DOMAIN WALL 

Now we consider the junction on Fig. 0] in which a ferromagnetic wire is inserted in between two superconductors. 
In the absence of a domain wall in the ferromagnetic wire (see Fig.0J-(a)) the junction is just a S/F/S junction. In 
the presence of a domain wall (see Fig.0J-(b)) Cooper pair-like states arising from crossed Andreev reflection can be 
transmitted through the junction. As a consequence the conductance is larger in the presence of a magnetic domain 
wall. We consider two limiting cases: 

(i) Transport is dominated by inelastic scattering in the ferromagnetic domains. Because of inelastic scattering the 
distribution functions in the ferromagnetic domains relax to the Fermi distribution. This case is discussed in 
sections HVBl ITvTl and HVDl 

(ii) Transport through the domain wall is phase-coherent and there is a Josephson current circulating between the 
two ferromagnetic electrodes. This case is discussed in section fVl 



A. The different time scales 



Similarly to Ref. 48 we notice that three time scales are involved in out-of-equilibrium transport through a ferro- 
magnetic domain wall: 

(i) The transport dwell time t c i being the time taken by an electron to travel through one of the magnetic domains. 

(ii) The energy relaxation time te- Because of inelastic scattering the distribution function in the out-of-equilibrium 
conductor relaxes to the Fermi distribution. This relaxation takes place on a time scale te- 
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(iii) The spin-flip time r s f being the time above which spin- flip scattering is relevant. 

We suppose in this section that t f <C t& <C r s r. The distribution function in the intermediate magnetic domains 
is thus well approximated by a Fermi distribution. The chemical potential of spin-up electrons is different from the 
chemical potential of spin-down electrons. 



B. Perturbative transport formula 

In this section we discuss the perturbative transport formula of the S/DW/N junction on Fig. EI The full transport 
formula to order i 4 is evaluated in Appendix The expression of I a ]l contains two kinds of terms: the terms (|B1|) - 
(|B3|I describe processes taking place locally at the interfaces between the superconductor and the ferromagnetic 
electrodes without propagation in the ferromagnetic electrodes. The terms l|B4j) - (|BT|) involve a propagation in the 
ferromagnetic electrodes. The two kinds of terms would contribute if electrodes (a, a') and (6, b') were single channel 
electrodes. We consider here that electrodes (a, a') and (b, b') are multichannel electrodes and we average the current 
over the microscopic phases. Once this averaging is done only the "local" terms survive in the transport formula given 
by 

-^XplfL [1 - PJ\ [n P (u ~ /V T ) - n F (u + M<a)] (38) 
^ 2 t 2 a t 2 p a p b {{g 2 a i/3 )) [1 + P a ] [1 + P b ] [n F {uj - fJ- a ,t) - n F {uj - MM )] (39) 
^ 2 tltl~p a p b ((fl p )) [1 + P a ] [1 - P b ] [n F (u - A*o, t ) " n F (u + p bd )] . (40) 

The term l|38|) corresponds to local Andreev reflection at the interface a - a. The term (I39|l corresponds to elastic 
cotunneling through the superconductor and the term (|40|l corresponds to crossed Andreev reflection. 

A similar calculation can be carried out at interface (a', a'). Once the average over the microscopic phase variables 



is carried out we find 

1% = -4 7 r 2 4,p a p'[l + P ][ riF ( a ;- Mo)T )-n F (cj-/x')] (41) 

+ SnHi (p a ) 2 {p'f [1 + P a ] 2 [n F (LU - Ma , T ) - n F {u; - //)] (42) 

+ SnHld^paPbdpl,^)) [1 + P a ] [1 + P b ] [n F (u - MM ) - n F (co - //)] , (43) 



where (I41|l and 142|) describe electron tunneling from the electrode (a, a') into the normal metal and l|43|l describes 
elastic cotunneling from electrode (b,b') to electrode (a, a'). 




C. Sequential tunneling through the S/DW/N junction 

In this section we discuss out-of-cquilibrium transport in a S/DW/N junction on the basis of the two-channel model 
shown on Fig. We suppose that a voltage V = is applied on the superconductor and a voltage V is applied on 
the normal metal. The spin-up and spin-down chemical potentials in the two magnetic domains (a, a') and (6, b') are 
determined in such a way that current is conserved. In general there are four unknown chemical potentials (/i Q ,Tj Mo,J.) 
Pb^l and Pb,i) that can be determined from four equations for current conservation. There exist two cases in which 
the 4x4 system of equations can be reduced to a 2 x 2 system of equations: 

(i) Half-metal ferromagnets where there is only one spin population in each of the ferromagnetic electrodes (a, a') 
and (b, b'). This case is treated in the main body of the article. 

(ii) The symmetric case where the two electrodes (a, a') and (b, b') have identical density of states and where t a = tp 
and t a i = tpi. This case is treated in Appendix IU1 

Let us consider half-metal ferromagnets: P a = 1, Pf, = — 1. The transport formula is found to be 

= ^m^P«AP*Ap'?^l,)) [1 - ^ 2 tl*Pa A p'] [1 - 2^t% Pbd p>] , (44) 

with 

V = t 2 jK'P^P'((flp)) t 1 - 2 ^>a,TP'] + 44*>*aP'<</a,/3>) [1 - ^%Pb,lP'] (45) 
+ t 2 a 4,(p') 2 [1 - 2^>a, T p'] [1 " ^Hl,p hA p'\ . 
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We note p_/v a typical value of the density of states either in the superconductor or in the ferromagnetic and normal 
metal electrodes. We first suppose that t^p^ and tap^ are small compared to t a i and tpi. The transport formula 
takes the same form as in the case where the ferromagnetic electrodes (a, a') and (6, b') are in equilibrium: 

^ = 16KXtlp aAPb , l {{fle)). (46) 
In the other limiting case where t^p^ and t%pjy are large compared to t a r and tp we find 



2 j-2 

107T p -5 5 ■ (4') 



itot 1C 2 / t a' t l3'PaAPb,i 

—— — lbn p -7. ~ 



We note g a = ld>-K 2 t 2 a , p a ,]p' and g b = l&ir 2 t 2 a , Pb,ip' the conductances associated to the interfaces (a', a') and (&',/?'). 
The total conductance is given by l/Gtot = l/<7a + 1/.96- The two interfaces are thus in series which is because 
transport is mediated by crossed Andreev reflection: a spin-up electron from the normal metal is transfered at site 
a', travels to site a, is reflected as a spin-down hole at site b. The spin-down hole travels to site b' and is transfered 
in the normal metal at site j3' . As a consequence of this transport process the two interfaces (a', a') and (b',/3 r ) are 
coupled in series. 

D. Sequential tunneling through the S/DW/S junction 

We consider the same model as in the preceding section but now the electrode on the right is superconducting (see 
Fig. EJ). We show that the properties of the S/DW/S junction are similar to the properties of the S/DW/N junction. 
We suppose that a voltage V is applied on the left electrode and a voltage V is applied on the right electrode. We 
consider a regime in which inelastic scattering in the ferromagnetic electrodes is strong enough so that the transport 
dwell time is much larger than the energy relaxation time (see section IIV All . Moreover we suppose that inelastic 
scattering is strong enough so that there is no Josephson effect. 

We consider that the ferromagnetic electrodes (a, a') and (b,b r ) are half-metal ferromagnets: P a = 1, P b = — 1. The 
total current is given by 

^tot 1e 2 ^a^^a'^P'({fa,fi))((fa',P')) s AO \ 

If we note g — 16n 2 t 2 i t'i((f^ /3 )) and g' — lQ^t^t 2 ,^/ 2 , /3 ,)) the conductances associated to crossed Andreev reflection 
at the contacts with each of the two superconductors we see that the total conductance is such that 

1 1 1 (49) 



Gtot g g' 

which is the expected result since electron pairs travel in series through the two superconductors. 

V. JOSEPHSON EFFECT IN THE S/DW/S JOSEPHSON JUNCTION 

A. A two-channel model 

Now we consider the S/DW/S junction on Fig. in which the right electrode is superconducting. We suppose 
that te ~ t<2 <C r s f or te < t& <C t s j so that the two superconductors can be coupled coherently through the two 
ferromagnetic channels. We look for the condition under which a Josephson current can circulate across the junction. 
We suppose in this section that the two electrodes (a, a') and (b, b') are half-metal ferromagnets with antiparallel 
spin orientations. The case of a partial spin polarization and non collinear spin orientations will be discussed in 
section IVBI 

The Nambu representation of the hopping matrix elements is given by 



t@, b = 



_£ ae i(V-X)/4 
t(je -i(<P+x)/4 

-tpe l ^ +x)/i 



(50) 
(51) 



11 




FIG. 6: Schematic representation of the Josephson junction containing a domain wall with a rotating magnetization. The x- 
and z-axis are shown on the figure. The y-axis is perpendicular to the figure. 



tb',/3' 



f Q ,e" i( >-x)/ 4 

-t^e^-x)/ 4 

t^e-^+x'/ 4 

-t p ,e i( - v+x ^ 4 



(52) 
(53) 



where (p is the difference between the superconducting phases in the right and left electrode and \ is the magnetic 
flux through the loop. We have the relations t QiQ = (4., a) , h,p — {tp,b) , i a >, a ' = (ia',a') and ipi^' = (tb',0') ■ The 
equilibrium current flowing from site a to site a is given by 



Ia,a = lj dun F {Lo)Tv [cr z [t a , a (g£ q - G R a ) - t a , a - G R a )] } . 

The Green's functions are 2x2 matrices since we do not discuss non collinear magnetizations for the moment. 
We deduce from 1541) that to order t the spin-up current through electrode (a, a') is given by 



(54) 



j(T) 



-2i-t a tpt a ,tp, simp 



A,ll A.22 
Ta.a' 9b.b' 



R.,11 R,22 
9a,a' 9b,b' 



Det[I-K A ] Bet[I-K R ] 
where K A ' R is the 4x4 matrix involved in the Dyson equation [I — K A ' R ]G A ' R — g A ' R : 



(55) 



l-K 1 ' 1 




-K 1 ; 1 

a ,a 


-K 2A 


-"■0,6 


1 — K 2,2 


^a'b 


-K 2 ' 2 

^b'b 


-K 1 ' 1 , 

a, a' 




i - kV , 

a' ,a' 


-if 2 ' 1 

^b'.a' 


^a.b' 


-K 2 ' 2 




i _ K ' 2 

1 ^b',b 





r C 1 ' 1 ' 




r a 1 ' 1 1 











U a,a' 




i.i 




L U a,b> - 








'2,1 
^ b,a 



= ^bfiU,' 



J-1,1 „1,1 „+„ 

■ L a,a9a,at CL, °- 



(56) 



The role of disorder can be included 



where we used the notation K^ a = t^ a g a<a t^ a g^, K b 

in a straightforward fashion. Since the spin-up and spin-down electrons of the Cooper pair propagate in different 

A X 1 A 2 2 

electrodes we should replace g a ' ; and g b b , ' by their average over disorder which decay exponentially with distance 
over a length scale equal to the elastic mean free patbS 4 -. We conclude that a Josephson current cannot be observed 
under usual experimental conditions since the size of the ferromagnetic region is usually much larger than the elastic 
mean free path. The opposite limit of small disorder is considered in Appendix iDl 



B. Multichannel effects 



1. Transport formula 

We consider in this section the multichannel S/DW/S junction on Fig. 0m which the ferromagnetic metal is multi- 
connected to the superconductors. The local magnetization makes an angle 9{z) with the z-axis. The situation with 
a uniform 9 corresponds to the multichannel 7r-junction. 
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The supercurrent is given by 



(57) 



where the Green's functions are 4x4 matrices. The supercurrent to order t 2 a t 2 p can be written as 



e r+ ° 
= -2i-t 2 a tpsm(p / ^ 2J /«m,a„//3 

" •'O 7, 1 _ _ 



(58) 



2, 2, A 1,1. A 
9g m .b k 9g n ,bi 

Det[/-i^] Det[/-if fl ] 



fc,Z,m,n 

2,2, fi 1.1. i?. 
9a m ,b k 9a n 'b, 



4,4, A 3, 3, A 
9a m ,b k 9a n M 



4,4,R 3,3,-R 
ia m ,b k 9a n ,bi 



Dct[I - K A ] Det[I-K R ] 



2,4, A 1,3, A 
9a rn ,b k 9a n ,bi 



2,4,R 1,Z,R 
9a rn ,b k 9a n ,bi 



T)ct[I - K A ] Bet[I-K R ] 



where o 

a a k 



is i-th Nambu component of the propagator connecting the two ends of the ferromagnetic metal at 
sites a k and b m . The first two terms in the r.h.s. of Eq. I|58|l correspond to a propagation without spin-flip in the 
ferromagnetic region whereas the last term correspond to a propagation with spin flip. We deduce from Eq. (|58|l the 
same conclusions as in the two-channel model. Namely an average Josephson current can circulate only if the size of 
the ferromagnetic region is smaller than the elastic mean free path in the ferromagnetic metal, a condition that is not 
usually realized in experiments. 



2. Limit of small disorder 

To obtain the supercurrent in the ballistic limit we replace the propagators by their expressions and sum over all 
channels. The propagators of a ferromagnetic metal with a rotating magnetization are not known in general. This is 
why we discuss here only the situation where the width of the domain wall is vanishingly small and the ferromagnets 
are half-metal ferromagnets. In this case there is no spin precession in the ferromagnetic region but there exist 
trajectories parallel to the interface that we can take into account in our calculation. The supercurrent is given by 

/ 2 \ 4 

/s = 8tt|z^ sin^f^J G(k F ), (59) 

with 

G(k F ) — — [° — / + °° — f° — [ +Q ° — f +0 ° dZm [ +OC dy,n f + °° — [ + °° — (60) 

■sin[k F R am:an ]sm[k F R l 3 k ,p l }cos[k F (R amybk - R anM )}, (61) 



2TrRa m ,a n 27Ti?/3 fe: /3 i 2irRa m ,b k 27rR a „,bi 



and with 



Ra m ,a n = \/ (zn - z m ) 2 + (y m - y n ) 2 (62) 

Rp k ,f3i = V \ z k - zi) 2 + (y k - yi) 2 (63) 

Ra m ,b k = VR 2 + (z m - z k f + (y m - y k ) 2 (64) 

Ra n< b, = VR 2 + ^n-z l ) 2 + {y n -yi) 2 , (65) 

where R is the longitudinal dimension of the junction (see Fig. [BJ. We have shown on Fig.[7|the variation of G(k F ) 
as a function of k F . We see that strong finite size effects are present but still we can make a comparison between 
(i) a calculation in which all trajectories are taken into account and (ii) a calculation in which only the trajectories 
perpendicular to the interface are taken into account. We see that for small values of k F (typically k F smaller than 
1/ao, where ag is the lattice parameter) the summation (i) is larger than (ii) whereas the opposite is true for larger 
values of hp. This shows that trajectories parallel to the interface play a relevant role in the determination of the 
supercurrent. 



VI. CONCLUSIONS 



To conclude we have presented a detailed investigation of several mechanisms involved in transport across several 
junctions involving ferromagnetic domain walls (S/DW, S/DW/N and S/DW/S junctions) The role of non collinear 
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FIG. 7: Variation of log 10 [G(kp)] as a function of kp. kp is measured in units of 1/ao- We restrict the summation in HilH to 
—L z /2 < z n ,zi < 0, < z m ,Zk < L z /2, —L y /2 < y n ,yi,ym,Vk < L y /2 and we use L — L y — L z . The open symbols correspond 
to all trajectories in IRjUl . The filled symbols correspond only to the trajectories perpendicular to the interface. 



magnetization was studied for the S/DW junctions. Using a spin ® Nambu (8 Keldysh formalism we have derived the 
form of lowest order transport formula valid for an arbitrary profile of magnetization. We find that the conductance 
is a scaling function of £o/D, where £o is the zero-energy BCS correlation length and D is the width of the domain 
wall. Because of the proximity effect an exchange field can be induced in the superconductor. Neglecting the spatial 
variation of the exchange field, we have derived the transport formula and shown that there was no spin precession 
around the axis of the exchange field. We discussed the transport formula of the S/DW/N junction. We have shown 
that to lowest order only the processes taking place locally at each interface played a role. These processes are: elastic 
cotunncling through the superconductor, crossed Andreev reflection, electron tunneling from the ferromagnet to the 
normal metal and elastic cotunncling through the normal metal. We described the transport of Cooper pairs across 
the S/DW/N and S/DW/S junctions in a regime where transport is dominated by inelastic scattering but spin is 
conserved. With these assumptions the local distribution function within the domain wall is a Fermi distribution with 
a different spin-up and spin-down chemical potential. This model provides a detailed description of the sequential 
tunneling of Cooper pairs across the S/DW/N and S/DW/S junctions. We described the Josephson effect in a 
S/DW/S junction. Diffusion is usually strong in a ferromagnet and disorder is thus expected to play a relevant 
role. In particular the Josephson current decays exponentially with the longitudinal dimension of the junction. The 
characteristic length is equal to the mean free path in the ferromagnetic metal. This means that a Josephson current 
cannot be observed in usual conditions. Nevertheless there can exist a finite current due to crossed Andreev reflection 
associated to elastic cotunneling in the ferromagnetic region. 
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APPENDIX A: SPIN PRECESSION IN THE METALLIC CASE 



In this appendix we consider the junction on Fig. [S] in which two ferromagnetic electrodes with non collinear 
magnetizations are connected to a normal metal£2*S. Our goal is to provide a comparison with the superconducting 
case presented in section All CI We suppose that a magnetic field h is applied on the normal metal and that the only 



effect of the magnetic field is to generate Zeeman splitting, 
elastic cotunneling takes the form 



The crossed conductance G a .b = 9I a /dVb associated to 



G a ,b — 8tt tatpPapb 



P a Pb sin 9 a sin 6^ cos I [&F,f — ^f.[ 



{1 + P a P b COS & a COs6 b 

Ra,p + eVb 



Spin precession can have two origins: (i) the term [fej?. 
a mismatch in the Fermi wave- vectors; (ii) the term 



l 

"-F.T 



«F4 
l_ 

«F, J. 



(Al) 



R a ,0 describes oscillations of the conductance due to 
R a .p describes oscillations in the conductance due 



1 


i 
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Ferro- 6 



a 



a 




Ferro- 6 b 



FIG. 8: The device considered in appendix^ Electrode ending at site "a" is a ferromagnet with a magnetization pointing in 
the direction 6 a - Electrode ending at site "b" is a ferromagnet with a magnetization pointing in the direction Oh- 
io a mismatch in the Fermi velocities. 

APPENDIX B: PERTURB ATIVE TRANSPORT FORMULA OF THE TWO CHANNEL S/DW/N 

JUNCTION 

In this appendix we provide a derivation of the transport formula of the S/DW/N model represented on Fig. (SJ 



The current through each link of the network on Fig. is given by the transport formula The spin-up current 
through the link a - a is found to be 



where and p a ,i are the spin-up and spin-down chemical potentials in electrode (a, a'), pb,l and pb,i are the 
spin- up and spin-down chemical potentials in electrode (b,b'), and p! is the chemical potential in the normal metal. 
After phase averaging we obtain three contributions to the transport formula: local Andreev reflection given by i|38fl . 
elastic cotunneling through the superconductor given by l|39|) and crossed Andreev reflection given by (|40|) . 



1. Transport at interface (a, a) 




~^ 2 t 4 a P 2 afL [1 - Pa] - Ma,t) - + Ma,|)] 

A^t^tppapbgi^ [1 + Pa] [1 + Pb] [n F (u - fi a ,t) -n F (u> - p b ,i)] 
AnHlfypaPbfl p [1 + P a ] [1 - Pb] [n F (uj - ti a>T ) - n F (uj + p b ,l)] 

47r 2 4*a'Pl,a'Pi,a'/loc i n F^ ~ /Vt) - n F {u + fl ad )] 

4:wH a t t a dp'pl a/ plyg ai 0g a ' t 0> [n F (uj - // ait ) - n F (uj - p b ,t)] 
4iT 2 t a t t a d 0/ pl al pl b ,f a ^f 0ta [n F {uj - p>a,i) - n F (u> + p b ,l)] 
4Trt a t0t a itp>Im gi*,gl$, g a ,0g a ',0'n F (uj - p'), 



(Bl) 
(B2) 
(B3) 
(B4) 
(B5) 
(B6) 



(B7) 



2. Transport at interface (a 1 , a') 



The same calculation can be carried out at interface (a', a'). The transport formula is found to be 




^ 2 tl4'P»Pb9a',0'9a',0' I 1 + P A [1 + P b] n F {uj - Mm ) 

8ir 4 tl,,t ,p a pbP 2 a >,(3> [1 + P a ] [1 + Pb] n F (uj - p'} 
A-kH^Pcp' [1 + P a ] [n F (uj - Pa,i) - n F (u - //)] 
8nH 4 a , {p a f (p'f [1 + P a ] 2 [u f (lu - /i a ,t) -n F (u- m')] 



(B8) 



+ 
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2mtitl,p'g (gl A a ) 2 n F {u - J) 

A 1A 
9a',p'9b,b' 



^■K 2 t a t t a it0ig a fip\ a ,Im 



+ 2mt a tpt a d &l g a ^p\ hl g£ >0 , g\ R a , + in p' 'g\ A a 



n F (uj - (j, a ,t) 

n F (u> - jU 6)T ) 



+ UlrtatptadpgafiPa!,? 9 a ,a' 9b,b' ~ in + P b)g a ,a> n F (uJ - fj,') 

- ATr 2 t 2 a ,tl,p a p b Re (5a',/3') 2 [1 + Pa] [1 + Pb] n F (w — p a ,t)- 

After averaging over the phase variables we obtain the transport formula given by Eqs. I|41|l - l|43[) that contains only 
two processes: tunneling from site a' to site a' and elastic cotunneling from site b' to site a' . 

APPENDIX C: TRANSPORT FORMULA OF THE SYMMETRIC TWO-CHANNEL S/DW/N AND 

S/DW/S JUNCTIONS 

In this appendix we consider S/DW/N and S/DW/S junctions with two symmetric channels. With this model we 
confirm the results obtained in the main body of the article for the asymmetric junction with half-metal ferromagnets. 
We suppose that the two channels have an identical density of states: p a = Pb, that the tunnel matrix elements are 
identical in the two channels: t = t a = tp, t' = t a > = tp> , and that the two channels have an opposite spin polarization: 
P a = P and Pb = —P. Then there exists a simple symmetry relation between the chemical potentials in the two 
ferromagnetic electrodes: p a ^ = pb.i and = pb,p 



1. The S/DW/N junction 

In the limiting case t <§; t 1 we have ~ V. The transport formula is identical to the case where electrodes (a, a') 
and (6, b') are in equilibrium: 



^ = 32^(1 - P 2 )/ 2 C + 32. 2 ty (1 + P 2 )«/^)). 



(CI) 



In the case of half- metal ferromagnets (P = 1) only the term corresponding to crossed Andreev reflection is non-zero 
and Eq. (|UT|> is equivalent to Eq. |gBJ|. 

In the limiting case pjvi' <C pNt <C 1 and t' <C pjvi 2 the current is the sum of a contribution due to local Andreev 
reflection and a contribution due to crossed Andreev reflection: J t ot = -Tar + ^car, with 



^AR 

'V 7 ' 

ICAR 



32^(trp P '(i-p 2 ).fL 


(<o> +<«£») 






/, 2 „c((^)) + ((/^)) 2 


+ a+p 2 )((f 2 , )) 





16^') 2 ,V«0} 


~(i + p 2 ) 


{((%,(,)) + <«/*») 


+ 2P 2 (fL- «%,?)))' 




'fLO&j)) + «/^» 2 ~ 


+ (i + P 2 )((fi f j)) 


'fL + ((9lJ)_ 





(C2) 
(C3) 



V 

In the case of half- metal ferromagnets (|C2j) and l|C3|l are equivalent to l|47|l. 

2. The S/DW/S junction 

In the case of the S/DW/S junction the total current is the sum of the local Andreev reflection and crossed Andreev 
reflection terms: 



1287T 



Iar 

V'-V V 

+ [tsai.p)) + {trail* [(i + p 2 )fL + (i - p 2 wi>, >))] } 



-^ 4 (i')V/i 2 oc(l - p 2 ) { [t\(gl, )) + {t')\{gl>M [(i - P 2 )fL + (i + P 2 Wl>M (C4) 
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^ry = 1 ^* 4 (* / )V<0)(i + P 2 ){^<fe)) + (i , ) 4 (<^^))] [{i-P 2 )fL + {± + P 2 Wl-,M (06) 

+ [* 4 «0> + (i') 4 «^,^»] [(1 + P 2 )fL + (1 - ^)<(/a' >/3 '))] } , 

with 

V = it s p 2 {(1 - P 2 ) [/^(s^)} + {{ftp)) 2 ] + (1 + ^P 2 ) [/ 1 2 c ((/ 2 , /3 » + (O>(<0Lj>>] } (06) 

+ 4(<') 8 P 2 {(1 ~ P 2 ) [/£ C «^,/J'» + ((Z 2 ',^)) 2 ] + (1 + P 2 ) [/l 2 oc((/ 2 ',^)> 2 + ((fa>,0>))((9t>,0<))] } 

+ u\t'f~ P 2 {(i - p 2 ) [fL{{gl, P )) + fL{{gl', >)) + 2«/ 2 /3 »«/ 2 ^»] 

+ (i + p 2 ) [fLdflp)) + fL{{&, >)) + ((fl, P ))((gl, P >)) + ((glM } ■ 

If the contacts with the two superconductors are identical we have t = t', f a ,p = fa',p' and g a> p — g a ' ,[3' from what 
we deduce 

IAR a 2 + 4 2 fl ( i u'2 



4^*V/£c(l-n (07) 



V'-V 

lr - vl = 47 r 2 tV<(^))(l + P 2 ), IC'M 



v- v 

where we used the notation /? = p/2 for the spin-up or spin-down density of state in the ferromagnetic electrodes. 
Eqs. (|C7|I and (|C8J| in the limit P = 1 are in agreement with Eq. 148(1 in the limit of a symmetric contact. In the 
symmetric case the conductance is thus equal to the conductance associated to a single superconductor divided by 
two, in agreement with Eq. (|49|l . 



APPENDIX D: JOSEPHSON EFFECT IN A TWO-CHANNEL BALLISTIC S/DW/S JUNCTION 

In this appendix we describe the Josephson effect with a ballistic propagation in the ferromagnetic electrodes. In 
the limit of a long junction R^> a the matrix / — K A ' R given by (|56[) is block-diagonal because the Andreev bound 
states do not couple the two superconductors. There exist two bound states associated to the interfaces (a, a) and 
(/?, b) and two bound states associated to the interfaces (a\ a') and (/?', b'). The secular equation for the bound states 
living at the interfaces (a, a) and (/?, b) takes the form 

1 + in 2 p F p N (t 2 a + t 2 ) ^ + (7r 2 p FPN ) 2 tlt f 2 ^ 2 ~ "° = 0, (Dl) 

where P n and P f are the density of states in the superconductor and in the half- ferromagnetic electrodes, and where 
we used the notation z — sin[kpR}/ (k F R). In the case t — t a — tp and in the tunnel limit irt 2 pNPF <C 1 the solution 
of I|D1|) takes the form 

uo 2 = A 2 [l + (Tr 2 t 2 PNPF ) 2 (l ± z) 2 ] . (D2) 

The supercurrent is easily deduced from 155|l : 

167r 9 e 4 2 2 2 a a Q a a . . . . r . , . 

(a o fc|0 2 (a o fcF) ft ^Ra,/3 2irK a >,p> ^R a ,a' ^Rb.v 

exp|— ^ ^ a ' a jl ^ h ' b ^ |/(z, z') cos a sin 9?, (D4) 



where a is defined by 



a = 4 (i? a ,a' - Rb,b> ) + 4" ( jR «.«' + ) ' ( D5 ) 



and where /(z, z') is a geometrical prefactor of order unity. 
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